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theoretical  and  experimental  determination  of  mechanical 

PROPERTIES  OF  SUPERCONDUCTING  COMPOSITE  WIRE*  \ 
W.  H.  Gray  and  C.  T.  Sunt 


INTRODUCTION 

The  object  of  tnis  research  is  to  characterise  the  mechanical 
properties  of  a  composite  superconducting  (NbTi/Cu)  wire  in  terms  of 
the  mechanical  properties  of  each  constituent  material.  In  1973,  the 
Cryogenics  Division  of  the  National  Bureau  of  Standards  (NBS)  published 
an  interim  report1  in  which  a  preliminary  investigation  of  the  mechan¬ 
ical  properties  of  a  solenoid  coil  composite  was  made.  The  coil  inves¬ 
tigated  consisted  of  epoxy,  fiberglass,  and  composite  superconducting 
wire.  Both  theoretical  and  experimental  elastic  constants  were  tabulated 


*  This  research  was  sponsored  in  part  by  the  Engineering  Researcn 
Iowa  State  Vniveralty,  toes,  IA  50011,  and  an  part 
by  the  Energy  Research  and  Development  Administration  under 
contract  with  Union  Carbide  Corporation. 

t  Department  of  Engineering  Sciences  and  Mechanics  and  Engineering 
Research  Institute,  Iowa  State  University,  Ames,  IA  50011. 
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for  a  typical  piece  of  coil  jput  out  of  a  small  solenoid.  Our  report 
differs  from  this  NBS  work  in  that  we  consider  only  the  mechanical 
properties  of  nn  individual  composite  superconducting  wire. 

The  theoretical  predictions  and  the  experimental  procedures  to 
determine  the  effective  elastic  constants  of  the  composite  wire  are 
described  in  the  next  two  sections  of  this  report. 

THEORETICAL  INVESTIGATION 

Most  of  the  analytical  work  for  predicting  the  mechanical  and 
thermal  properties  of  fiber-reinforced  composites  in  terms  of  volu¬ 
metric  composition,  geometrical  arrangement  of  the  fibers,  and  con¬ 
stituent  material  properties  was  done  before  1970.  There  are  five 
approaches  to  predict  the  micromechanical  behavior  of  fiber-reinforced 
composites.*  The  essential  characteristic  of  each  is  described  below. 

1.1  SELF-CONSISTENT  MODEL  METHODS 

3  A 

This  method  was  originally  proposed  by  Hershey  and  Kroner  for 
crystal  aggregates,  and  was  first  employed  by  Hill^  to  derive  expres¬ 
sions  for  elastic  constants.  Hill  modeled  the  composite  as  a  single 
fiber  embedded  in  an  unbounded  macroscopically  homogeneous  medium, 
subjected  to  a  uniform  loading  at  infinity.  This  uniform  loading 
produces  a  uniform  strain  field  in  the  filament  which  is  then  used 
to  estimate  the  elastic  constants.  A  similar  model  proposed  by  Frohlich 
and  Sack6  for  predicting  the  viscosity  of  a  Newtonian  fluid  containing 
a  dispersion  of  equal  elastic  spheres  consists  of  three  concentric 
cylinders,'  the  outer  one  being  unbounded.  The  innermost  cylinder  is 
assumed  to  Have  the  elastic  properties  of  the  filaments;  the  middle 
one  has  the  properties  cf  the  matrix;  and  the  outermost  has  the  properties 
of  the  composited  The  solid  is  subjected  to  homogeneous  stresses  at 
infinity.  The  resulting  elastic  fields  are  determined,  and  then  are 


*  No  attempt  is  made  here  to  give  a  comprehensive  literature  survey 
regarding  this  subject.  More  references  can  be  found  in  Ref  •  2. 
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employed  to  predict  the  elastic  constants  of  the  composite.  Applications 
of  the  self-consistent  model  methods  can  be  found,  for  example,  in 
Refs.  7-9. 

1.2  VARIATIONAL  METHODS 

In  this  method,  the  energy  theorems  of  classical  elasticity  are 
used  to  obtain  bounds  on  the  mechanical  and  physical  properties  of 
filamentary  composites.  The  minimum  complementary  energy  theorem 
yields  a  lower  bound,  while  the  minimum  potential  energy  theorem  yields 
the  upper  bound.  Using  this  approach,  bounds  for  the  elastic  and 
thermal  properties  of  composites  have  been  obtained  by  many  investiga- 


1. 3  EXACT  METHODS 

By  assuming  that  the  fibers  are  arranged  in  a  doubly  periodic 
rectangular  array,  a  fundamental  or  repeating  element  can  be  estab¬ 
lished.  The  resulting  elasticity  problem  can  then  be  solved  either 
by  introducing  a  stress  function  using  a  series  development,  or  by 
numerical  techniques  such  as  finite  difference  or  finite  element 
methods.  Once  the  problem  is  solved  elastically,  the  resulting  elastic 
fields  can  be  averaged  to  get  expressions  for  the  desired  elastic  con¬ 
stants.  Typical  applications  of  this  method  can  be  found  in  Refs. 

13-18. 

1.4  MECHANICS  OF  MATERIALS  METHOD 

By  making  simplifying  assumptions  regarding  the  mechanical  or 
thermal  behavior  of  a  composite  material,  the  mechanics-of-materials 
expressions  for  the  equivalent  elastic  or  thermal  constants  of  unidi- 
rectionally  reinforced  fibrous  composite  materials  can  be  derived.  For 
example,  to  determine  the  longitudinal  Young's  modulus,  one  assumes 
that  the  longitudinal  strains  in  both  the  matrix  and  the  fiber  are  the 
same;  in  order  to  determine  the  transverse  Young’s  modulus ,  one  assumes 
that  transverse  stresses  in  both  materials  are  the  same.  This  approach 
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usually  is  referred  to  as  the  "rule  of  mixtures."  The  "rule  of  mixtures 
expressions  for  elastic  moduli  and  thermal  conductivities  can  he  found  in 
Refs.  19-21. 

1.5  THE  HALPIN-TSAI  EQUATIONS 

For  designers,  it  is  often  necessary  to  have  simple  and  rapid  com¬ 
putational  procedures  for  estimating  the  macromechanical  properties  of 
a  fibrous  composite.  Such  empirical  formulas  have  been  developed  by 
Halpin  and  Tsai22  based  upon  modifications  of  the  results  discussed 
under  approaches  1.1  and  1.3.  By  estimating  the  value  of  a  factor 
which  depends  on  the  geometry  of  the  inclusions,  spacing  geometry,  and 
loading  conditions,  the  composite  elastic  moduli  can  be  approximated. 
Reliable  estimates  for  this  factor  can  be  obtained  by  comparing  the 
Halpin-Tsai  equation  with  the  numerical  micromechanics  solutions,  -f 
used  appropriately,  the  Halpin-Tsai  equation  can  yield  very  reliable 
results  without  elaborate  calculations. 

All  the  above  methods  make  the  following  three  basic  assumptions: 

(1)  each  constituent  material  behaves  linearly  elastically,  (2)  the  fibers 
are  straight  (without  twist),  and  (3)  there  are  no  residual  stresses. 

For  this  investigation,  we  use  several  of  the  available  theoretical  equa¬ 
tions  to  predict  the  effective  elastic  mechanical  properties  of  supercon¬ 
ducting  composite  wire.  These  equations  are  listed  in  Appendix  I. 

In  general,  a  superconducting  composite  wire  may  be  twisted  to 
minimize  ac  power  losses  in  a  superconducting  magnet.  Twisting  of  the 
wire  violates  an  assumption  implicit  in  the  derivation  of  all  the 
equations  presented  in  Appendix  I.  However,  we  believe  that  this  effect 
is  small  (see  Appendix  II),  and  for  engineering  purposes  can  be  neglected. 
Other  effects,  such  as  inelastic  behavior  of  the  wire  at  higher  loading 
levels,  and  residual  stresses  in  the  wire  introduced  during  fabrication, 
may  influence  the  results  presented  in  this  paper,  and  should  be  ana¬ 
lyzed  more  thoroughly. 
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2.1  EXPERIMENTAL  DETERMINATION  OF  En  AND  V12 

^11  aru*  V12  can  determined  from  a  simple  tension  test  (see 
Fig,  1),  The  direction  of  loading  is  parallel  to  the  fibers  of  the 
conductor.  The  longitudinal  Young's  modulus  is  ^.‘ermined  from  a  0^ 
vs  diagram  where  is  equal  to  the  applied  load  divided  by  the 
cross  sectional  area  of  the  specimen,  and  is  the  strain  along  the 
fiber  direction  of  the  conductor*  Figure  2  represents  an  experimentally 
determined  plot  of  this  diagram  for  KRY0-210  superconductor  showing  a 
value  for  E^  of  119  GPa.  A  comparison  of  the  theoretical  prediction 
and  the  experimental  data  is  shewn  in  Fig,  3  in  this  graph,  as  well  as 
the  four  normalized  comparison  graphs  which  follow;  the  legend  refers 
to  the  theoretical  equations  presented  in  Appendix  I,  Both  equations 
predict  the  same  behavior,  which  deviates  approximately  3%  from  the 
experimental  data. 

The  major  Poisson's  ratio  is  determined  from  the  slope  of  the 
va  diagram  during  the  same  experiment,  where  is  the  strain  in 
either  transverse  direction.  The  experimentally  determined  value  was 
0.347  (sec  Fig.  4).  The  normalized  plot  (see  Fig.  5)  showing  the 
comparison  between  theoretical  prediction  and  experimental  value  again 
demonstrates  little  difference  between  theories  with  the  experimental 
data  differing  from  the  predictions  by  about  2%. 

2.2  EXPERIMENTAL  DETERMINATION  OF  E22  AND 

E22  and  are  determined  from  a  test  similar  to  that  described  in 
2.1  (see  Fig.  6).  The  direction  of  the  applied  load  is  normal  to  the 
fiber  axis.  The  transverse  Young's  modulus  is  determined  from  the  cr2 
vs  e2  diagram  where  02  is  the  stress,  and  e2  is  the  strain  in  the  direc¬ 
tion  of  the  applied  force.  Figure  7  represents  an  experimentally  deter¬ 
mined  plot  of  this  diagram  for  KRYO-210  superconductor  showing  a  value 
for  E22  of  122  GPa.  The  experimental  data  are  compared  to  the  theoretical 
predictions  for  E22  in  Figure  8.  An  error  of  approximately  5%  is  observed. 

Analogously  the  minor  Poisson's  ratio,  and  therefore  G23,  is  deter¬ 
mined  from  an  vs  e2  diagram.  The  experimental  value  for  G23  is  43.1 
GPa  (see  Fig.  9).  Figure  10  compares  this  data  point  with  the 
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.Helical  predictions.  An  error  of  approsi.ar.iy  2,  is  observed. 

2  3  EXPERIMENTAL  DETERMINATION  OF  THE  LONGITUDINAL  SHEAR  MODBUS  C,, 

The  strap  lest  way  to  derer.lne  the  longitudinal  ^ 
is  to  ose  a  rensite  specimen  -h  the  fibern  ^ 

ro  the  geometrical  axis  of  the  spec»e».  ^ 

a  composite  sopercondoctor,  however,  would  be  verj  dittl 
to  fabricate.  An  alternate  method  to  evaluate  C,2  experlmen 

'““tmTiwo-diraensionel  anisotropic  stress-strain  relation,  we 
have 


0 
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?16 

*26 
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in  the  xz  plane  and  make  angles  o  ana  o  -r  ^ 

tively;  Y,i  is  the  shear  strain  of  x'  and  z  axes,  a  6’  0  +  */2 

,«  tie  corresponding  normal  and  shear  stresses  respectively  (see 
T®  .  fci  4  Eq  (2)  represents  the  compliance  it  .trix 

Fieure  11).  The  matrix  [SJ  xn  fcq.  u; 

ini  tolposite  .terUl  in  xV  directions.  In  a  simple  tension  «- 

with  the  applied  load  and  the  fibers  oriented  along  the  X  ”  ’  ® 

and  te  4  can  be  measured 

data  obtained  from  strain  gage  ros 
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where  is  equal  to  the  applied  force  divided  by  the  cross  sectional 
area  of  the  specimen*  The  elastic  compliances  S^,  S 12>  S^,  S  % 

and  ^66  are  related  to  En>  E22*  vio*  and  Gi?  by  the  following  relations.2 


-  cos40 

*11  E  +  *■ 

XI  En 


Mr* 


sin20ccs2O  + 

22 


(4) 


5J2  -fr^-  +  ■ET - sin2Ocos20  -  (sin  *0  +  cos 

\bll  E22  G12/  hi 


cos  0) 


(5) 


(6) 


®22  -  +■  2(<r-  -  i^sinWo  + 

22  En  'G12  E11  /  E22 

'  (ft  +  ft  '  ft )  *  (ft  *  ft  -  4)  5l“3e“5°  "> 

?26  '  (ft  *  ft  ‘  ft)  81,Ac0sS  +  (ft  '  ft  "ft)  Sl“0'o53°  <8) 

2  v 

S66  “  (£  +  +  )  sin20cos20  +  (sin40  +  cos40)  (9) 


Since  E^,  E22,  and  v12  are  determined  from  tests  (2.1)  and  (2.2),  and 

e6*  e9  +  v/2*  Y0*  a0’  Q0  +  tt/2*  and  T0  are  obtaiRed  either  from  direct 
measurement  by  strain  gages  or  from  Eq.  (3),  the  only  unknown  in  Eq.  (2) 

is  C12.  Thus  G^2  can  be  computed  from  any  one  of  the  three  equations  in 
Eq.  (2).  Its  value  is  found  to  be  44.8  GPa  which  is  within  5%  of  the 
theoretical  predictions  (see  Fig.  12). 


CONCLUSIONS 

The  goal  of  this  experiment  was  to  determine  which  theory  of  com¬ 
posites  best  predicted  the  elastic  mechanical  behavior  of  a  supercon¬ 
ducting  (NbTi/Cu)  composite  wire.  Examination  of  each  elastic  mechanical 
property  reveals  that  all  theories  examined  are  capable  of  predicting 
experimental  data  to  within  5%. 
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Sln«  the  Poisson's  ratios  tor  both  the  fiber  and  the  matrix  tat.rlal 
were  very  close,  there  was  essentially  no  (less  than  1%)  difference  among 
all  the  theoretical  predictions  for  an.  Individual  neohanlcal  constant. 
Because  of  the  expense  and  difficulty  of  producing  elastic  constant 
within  0.1*  accuracy,  and  therefore,  conclusively  determining  which 
theory  Is  best,  no  further  experirmnts  were  performed. 

in  conclusion,  for  a  supercondu.ting  composite  wire,  hb.i/Cu,  a 
simple,  fast,  and  reliable  engineering  estimate  of  its  elastic  mechan  c. 
Ltv L  can  be  made  by  tming  the  "rule  of  mixtures."  It  is  unnecessary 
M  „se  one  of  the  more  rigorous  theories. 
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NOMENCLATURE 


Symbol  Definition 

E  -  Young1 s  modulus 
G  -  Shear  modulus 

V  -  Poisson's  ratio 

K  -  Bulk  modulus  in  plane  strain 

V  -  Volume  fraction 

f  -  Fiber  material  (NbTi) 
m  -  Matrix  material  (Cu) 

L  -  Composite  material  constants  along  the  fiber  direction 
T  -  Composite  material  constants  along  a  transverse  direction 


APPENDIX  I 


This  appendix  presents  the  various  theoretical  equations  which  were 
used  to  predict  the  elastic  mechanical  properties  for  the  superconducting 
composite  wire  studied  in  this  report.  Their  alphabetic  index  refers  to 
the  legend  on  the  particular  graph  where  a  comparison  with  experimental 
data  was  performed.  Also  included  in  this  appendix  is  the  reference 

in  which  these  equations  may  be  found. 

A  program  which  numerically  tabulates  all  of  these  equations  is 

listed  in  Appendix  III. 
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1.  Longitudinal  Youngs  modulus 

(a)  E  =  EV  +  EfV- 
v  7  L  mm  ff 

Source:  Reference  2 


<W  el  -  (y. ♦  e£ve)  » 

Em  (D1  ~  D3F1^  1  E£^P2  ~  D4F2) 

*  =  V°1  "  V  +  E£^°2  -  V 

1  +  Vr 

D,  -  1  -  V- 


D2  = 


+  V 


m 


D3  =  2Vf 


2  Vf 

°4  “  2  V  “ 

m 


V  V,E,  +  V,v  E 
_  m  f  £  £  m  m 

F1  +  V  E  V- 


Tff 


L-  V  £ 

m  f 


F2=-Fl 

m 


(10) 


(11) 


Source 


Reference  10 


2 


Ma-jor  Poissonts  ratio 


VfEfL1  -1-  VmEmVmL2  (12) 

(a)  VL  =  VfEfL3  + 

■  4  -  2Vf(l  -  v„2)vf  +  v.(l  *  v„)f„ 

t2  -  (l  - .  V,  - 

i3  -  2(1  -  v„2)vf  +  (l  *\)\ 

Source:  Reference  10 

(b)  VL  *  vmVm  +  VfVf 

Source*  Reference  2 


(13) 
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Source:  Reference  23 


4,  Transverse  shear  modulus 

r  2(1  -  v  ) 

(a)  Gt  (lower  bound)  «  cl  1  -  T  J  v  '  VfA4 

*  Vm 


£ 

where  is  obtained  by  solving  the  following 
equations 


Source:  Reference  10 


and  A?  is  obtained  by  solving 
4 


Source:  Reference  10 
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(24) 


(25) 


(26) 


appendix  II 


Whitney25  has  investigated  the  influence  of  twist  on  graphite 
in  an  epoxy  natrix.  He  derived  an  equation  for  the  reduction  of  the 
longitudinal  elastic  Young* s  *>dulus  for  graphite  as  a  function  o  he 
geometry  of  the  fibers.  This  equation  is  directly  applicable  in  esti¬ 
mating  how  twisting  affects  a  superconducting  wire. 

If  the  initial  and  reduced  moduli  are  Ej  and  Ep,  respect  ve  y, 

.  .25 

their  ratio  can  be  expressed  as. 


EI  1  +  4ir2N2R2  , 

«h«r.  »0  is  the  number  di  twists  per  centimeter  «nd  R  is  the  radius  •£ 
the  fiber  in  centimeters. 

For  the  superconducting  wire,  analyzed  in  this  report. 


N  *  .132  cm 
o 


-1 


„-3 


R  »  3.17  *  10  cm, 
which  yields 

1  _  _ 


2  2  7. 
1  +  AttVr 


1  +  6.87  x  10 


.-6 


=  1. 


Clearly,  the  effect  of  twist  for  this  superconducting  wire  can  be 
neglected. 


23 


\ 


APPENDIX  HI 


I 


26 


uuo 

Suu 
uu 
uuu 
uuo 
uuu 
uuo 
uuu 
uuu 
uuu 
uuu 
uuu 


uouuoouuoooouuu 

UOUOOOUOUOUUOOO 

uuuuuuouuuuuuuu 


uuu 

uuu 

uuu 

uuu 

uuu 

uuu 

uuu 

uuu 

uuu 

uuu 

uuu 

uuu 


V? 

CM 

CD 

u> 


oca  OT 
ocacoc 

0CQCQC 

ococoe 

OCOCQC 
dc  a  a 
OCOCQC 
DC  DC  DC 

ocacoc 

OCOCQC 

ocacoc 

OCOCQC 


ocacoc  acococ 

OCOCQC  OCOCQC 

ocacoc  ocacoc 


ocacoc  acococ 
acococ  acococ 
acococ  ocacoc 


OCOCQC 
OCOCQC 
OCOCQC 

OCOCQC  OCOCQC 

OCOCQC  OCOCQC 

acococ  acococ 

OCOCQC  OCOCQC 

acacacacaca: 

OCOCQC  OCOCQC 

ocacoc 

aocac 

ocacoc 

QCQcQCOCQCaCQCOeaCOCQCQC 

QCQCQCQCQCQCOCQCQCQCaCQC 

QcacQCQcacQCQcacQcacocac 


uuu 

lililiUi 

UUU 

UlliJUi 

uuu 

uuu 

tljuitil 

UUU 

UUU 

UUU 

UUU 

UUU 


UJUItU 
UUU  til 
UIUJUI 
UJUJUJ 
UIUIUI 
UlUJUi 
UJUJUJ 
UIUIUJ 
UIUIUI 


UJUIUIUJUJUIUJUIUUJUJUJUJ  UlUJUi 

'jj  l|ILil  / 1 1  111  I  ■  I  111  111  111  111  111  LlJ  1  i  i  1 U  111  LlI 

UUUUUUUUUUUUUUUU 


UUIUI 
UUJUI 
UILUUJ 
UUJUI 

uuu 
uuu 
uuu 
uuu 
uuu 
uuu 
u  u  u  uu 
uuuuu 

U  U  U  U  UJ 


X 

ac 

o 


U)_J 

tz 

3*OC 

£§ 

I  « 


c  n 

—  M 


?  »n  \ 


CD  ffi  CD  CD  CD  CD 

QJ  (Q  QJ  (])  fj)  0) 

ffiCOffl 

acorn 

CD  CD  CD 
comm 
moon 

CD  CD  CU 
CD  CD  CD 
CD  CD  CD 
CD  CO  CD 
CDCDCDCQOQCDCDCDCD 
CDCDCDQICDCDGQCDCD 

CO  CD  GD  CD  CD  CD  03  03  QD 


CD  CD  CD 
CD  CD  CD 
CD  CD  CD 
CD  CD  CD 
CD  CD  CD 
CD  CD  CQ 
CD  CD  CD 
CD  CD  CD 
CD  CD  CD 
CD  CO  CO 
CD  CD  CD 
CD  CD  CD 


fC  00  CD  CD  (Q 
CD  ^  ^  CD  CO 


Cl)  CD  CD  CD  CD 
CD  CD  QD  QD  QD 
QD  CD  CD  (D  CD 


CD 

CD  CD  CD 
GO  CD  CD 
QD  CO  CD 
CD  CO  CD 
CD  CD  ID 
CD  CD  CD 
CQ  CD  CQ 
CD  CO  00 
CD  CD  CO 
CD  CD  CD  CD 
CD  CD  CD  CD 
CD  CD  CD  CD 


ittaaaaaaffcc 
a<r<iaa<rctaaa:<x 
aaaaaaaaacccc 
a  a  a  era 

aaa  a  a 

a  a  a  a  a 

arcc  a  a 

aaa  a  a 

Saa  a  a 

a  a  a  a 

a  a  a  a  a: 

caa  a  a 

aaaaaaaaaaa 
aaaaaaaaaaa 
aaaaaaaaaaa 


aaaaaaa 

aaaaaaa 

aaaaaaa 

a 


a 

aaaaaaa 

aaaaaaa 

aaaaaaa 


lU.IL  U.U.U. 

.  u_  u.  u.u.u. 

_  U.  U_  ILU.U. 

lU.Il  U.  U.U. 

lU.il  uuu 

LU.U.  ILU.U. 

„UU  UUU 

-  U.  LL  U.U.U. 

.  U  U  MLU 

.UIlUU-U-1L.UUU.LlLlU_IJ_LLLlU.UU.U-LL 

»UUUUULlUUU_UUUUUUUUUUU 

.lLtLlLlLlLU.U.U.U.lLtLlLU.lLU.lLlLlLU.U. 


ou 

va 

-o. 

ID  — 
k  - 
t  w 
uu 
C  — 

a  a. 
i  o 


su 

CD  DC  -  — 

ou  u—  r 
jm  mu  'j 

a.  —  -'to  - 

UU 

O  M  U 
C  OLU  ft 
O  JD  -  u  J  — 
O^T  ~  O 

^L,  L. 

c  a 

-nyio 
CKl  — 

CCD  NON 

d  j*i  »  a  in 

acvuoo:  CD 
«Moav 
CD  *  CM - 
«->CD  O 

<Noaz  -» 

|U|>£CD 

wa  a  a 

u  oe  >•  u  x 

uum«uj 
c  »—  woe 

O  i.  -'ll.  0 
-OrC«£# 
CHIiDUfl 
LD  L  CD 

O 

>  *  1/1  3  — 

-'CH/)- 

ja  n  3 

o.  a  o  -u 
cm-  3  od  o 
F-  cn  cr— uj  — 
cl  *  ®  -  n  — 
j  auoiL 


noonoononnnoonoo 


27 


PROGRAM  FIBERC.V1A 

LAST  UPDATE’  ,  IAPR76 

■j 

AUTHOR:  U.  .<•  GRAV 

0AK°RIDGE  NATIONAL  LABORATORY 
OAK  RIDGE.  IN.  37630 

LANGUAGE:  DECSYSTEM-18.  FORTRAN- 18 

SUBROUTINES  REQUIRED: 

HINV  -  STANDARD  IBM  SCIENTIFIC  SUBROUTINE  PACKAGE. 

IHHICIT  REAL  «> 

REAL  NUF.NUli  ... 

DlfENSlON  UC6#6).RES<6)*L(b)#n<6>pHtoJ 

•  d'HnIIoN  ELEM(2)^>WUt.T<2)»GLGM(?>«GTGH(5>» 

‘  data'pi^SiIw^ 

DATA  IORED. IOURT/3.5/  ■  < 

c  BEGIN  PROGRAM  F1BERC.VIA 
C  OUERY  THE  USER  FOR  THE  MATERIAL  PROPERTIES 
„  $$m?*TFIBERC.VlA  DOCUttNTED  IN  0RHL/TO-S331V 

11  I  'INPUT  IN  THE  FOLLOWING  ORDER:*/ 

2  •  EFCFIBER  YOUNGS  MODULUS)  / 

3  #  EMC MATRIX  YOUNGS  MODULUS)  ^ 

4  *  NUF (FIBER  POISSONS  RATIO)  / 

\  : 

l  '  ReSIweI  A^ SP ACE%EL*I MI TS  the  INPUT  VARIABLES-^) 

0  reSduored.i8>ef.em.nuf.;;um.vf 
10  FORMAT  (5F) 

909  CONTINUE 

C  CALCULATE  THE  SHEAR  «DULII 
GF-EF/(2.*< l.+NUF)) 

GM-EfV(2.*U.^NUM))  COftrTfnM 
C  CALCULATE  THE  MATRIX  VOLUTE  FRACTION 
VM-l.-vf 

C  CALCULATE  THE  BULK  WDULM 
KF -F.F/C  3 .  *U . -2 .  *«UF  >  > 

KM»EM/<3.*(  1**"2.*NUM) ) 

C  CALCULATE  THE  FIBER-MATRIX  RATIOS. 

GMGF  *GM/GP 
GFGH-GF/CM 
EFEM-EF/EM 
EMEF-EIVEF 
Gf€M*GM/EM 
KMGM-KM/GM 
KFGF-KF/GF 
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KFGM-KF/G M 
C  ECHO  INPUT  DATA 

URITEU0URT*29>  EF#EM*NUF*NUM,VF,GF*GM*EFEM,KF#KM*GFGM 
29  FORMATC//*  INPUT  DATA*//* 

I*  EF  -  **T10*IPE11.4*/. 

2'  EM  •  **T10#  IPEU.4#/# 

3*  NUF  -  * • TIB# IPEI 1.4*/* 

4*  NUfl  -  * #T10#  tPEI  1 .4#/* 

5*  VF  •  *#T10* 1PEI1.4#/* 

6*  GF  •  **T10*  IPEU.4,/# 

2*  GI1  •  **T10* IPEl 1.4#/, 

8*  EF/Efl  -  **T10#  IPEl  1 .4*/* 

9*  iCF  •  *#T10*  IPEl 1.4#/# 

.*  KM  -  *#TI0#  1PEI1.4#/#*  GF/GM  •  **T10#  1PEU  .4) 


C 

C 

C 

c 

c 

c 

c 

c 


21 


20 

25 

C 

eei 


38 

24 

C 

c 

c 

c 


THE  EQUATION  NUMBERS  GIVEN  IN  THE  REST  OF  THIS  CODE  REFER 
TO  ORNL/TH-5331 

CALCULATE  THE  LONGITUDINAL  YOUNG'S  MODULUS 


EQUATION  10 

ELEHCl)*VF*EFEM*VM 
EQUATION  U 

IFCVM.NE.0.0)  GO  TO  21 
ALPHA* 1. 

GO  TO  20 
CONTINUE 

F  l  •  CNUM*VF*€FEMH«JF*VM)  /CNUF*VF*£FEM+NUF*VM) 

F2*NUr«Ft/NUM 

01*1. -NUF 

D2-<l.+vF>/VM+NUM 

D3*2.*NUF*NUF 

D4-2 .  *NUM*NUM*VF/Vt1 

ALPHA*Dl-D3*Fl4£FEM*CD2-D4*F2> 

ALPHA*ALPHA/(D l-D3+€FEtt*( D2-D4) ) 

ELEMC2) -ELEMC l)*ALPHA 
UR I TEC I0URT*25) 

FORMAT (//# *  LONGITUDINAL  MODULUS*/) 

OUTPUT  DATA  ON  DATA  FIlE 
DO  601  J*l#2 
ELEMC  J) -ELEMC  J)*€M 
DO  30  J*I#2 
JJ-J+9 

LRITEC  IOLRT#24>  JJ*ELEM(J> 

CONTINUE 

FORMAT ( IX* 'EQU*  *# 12* 4X. IPEl 1.4) 


CALCULATE  THE  MAJOR  POISSON  RATIO 


EQUATION  12 

>0.1*2.  *NUF*C  I .  -NUM*NUM>  *V^4NUM* ( 1 . 4NUM)  *VT1 
XL2*VF*1 1 . -NUF-? . *NUF*NUF) 

«-3*2 . *C 1 . -NUM*NUM) *VF+C 1 . +NUM) *VM 
XNULTA-VF*£FEM*XL 1 +VM*NUM*XL2 
WULTC  l )  *XNULTA/C  VF*€FEM*XL3«-VM*rt.2> 


C  EOUATION  13 

XNULTC2) -VF*NUF>VM*NUM 
URITEUOURT*  225) 

225  FORMATC//.*  MAJOR  POISSONS  RATIO*/) 
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■  OUTPUT  DATA  OH  DATA  FILE 
DO  31  J*l»2 

URITEUOIRT.24)  JJ»>NULT<J> 

It  CONTINUE 

CALCULATE  THE  TRANSVERSE  YOUNG'S  WCULUS 

5  iouatiohw^ >/CVWB€Wr> 

C  EOUAT^ON  15  J 

ETEMC2) *1  I  .♦2.*€TA*VF)/t  I  .■ETA«'^> 

l  EQUATION  .6  ,HUFJ) 

Sh-1  ./<2.*<  I  .-HUM)  >  ^  _ 

‘  E0WW.2..<E*F-I.» 

C  CHECK  FOR  TO  160 

|tEN( 4?  •  1 1- 2*.  »SQRT^v?^P I ) ♦ ( P I -SO*TANSTF ) /ALPHA 

GO  TO  161 

160  CONTINUE 
ETEM(4>  —  i. 

161  CONTINUE 

*  E^0^EFE^KUr)«2*VE^ 

ETEMC5)  -Vt1^mVF/EF-S0/(VF<€F'"V^trHl . ) 

ETEM(5)  •  I .  '(ETEMC5)  *€M> 

625  F^RMATl'^* TRANSVERSE  YOUNGS  MODULUS''* 

C  OUTPUT  DATA  ON  DATA  FILE 
C  EHG-Ar*Xt<EF,EH> 

ELU*AMN1CEF,EM) 

DO  B65  J»1#T- 

005  ETEMU^ETEMCJ)*^ 

DO  32  J*l#5 

IF(ETEM(J) .GT.EHG.ORjETEHCJ) .LT.ELU)  CO  TO  02 
URITE(I0LRT#24)  JJrETEMCJ) 

GO  TO  32 

H  ®TCS™  HOT  APPLICABLE') 

32  COHTINUE 


02  URITEt 

03  FORMAT 

32  COHTINI 

C  CALCULATE  T 

C  EQUATION  19 
C 


CALCULATE  THE  TRANSVERSE  SHEAR  IDDULUS 


UU  K.OO 


// 

I 

i 
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C  BEGIN  TO  SET  UP  THE  COEFFICIENT  MATRIX 


U(M) 

•l.  * 

U(2. 1) 

-0. 

0(3.1) 

-1.  1 

U(4. 1) 

•0. 

0(5.1) 

-1. 

U<6. 1) 

•0. 

0(3.2) 

•1.  , 

0(3.3) 

•1.  1 

UC4.3) 

-2.  i 

U15.3) 

-3. 

0(6.3) 

-2. 

0(3.4) 

■1. 

U(  1.5) 

-8. 

0(2.5) 

-8. 

0(3.5) 

—  1. 

0(4.5) 

-0. 

0(6.5) 

•0. 

U<  1.6) 

•0. 

0(2.6) 

•0. 

U(3.6) 

—  1. 

0(1.2) 

•I.AF 

UC4.2) 

—  (3.-4.*NUM)/<3.-2.*NUM> 

0(2.2) 

•0(4.2)^vF 

UC6.2) 

— l./(3.-2.*NUM) 

U(5.2) 

—3. *0(6.2) 

0(1.3) 

•VF*VF 

0(2.3) 

— 2.*un.3> 

U( 1*4) 

-vE 

UC4.4) 

•I./(1.-2.*NUM) 

U(2.4) 

•VF*U(4.4> 

0*5.4) 

•0(4.4) 

0(6.4) 

— U(4.4) 

U(5.5) 

•-GFGM 

0(4.6) 

«(3.-4.*NUF)^(3.-2.<WUF) 

UC6.6) 

•GFGM/(3.-2.*NUF) 

0(5.6) 

•-3. *0(6. 6) 

STORE  THE  COEFFICIENT  MATRIX  INTO  UINV 
DO  71  IH-1.6 
DO  71  Jrt-1.6 
UlNVl  IW. JH)  »U(  IH.  JH) 

FORMULATE  RIGHT  HAND  SIDE 
AUJ-l. 

AC2)-0. 

A<3)>0. 

A(4)-0. 

A<5) »8. 

A(6>-0. 

C  INVERT  THE  COEFFICIENT  MATRIX.  SUBROUTINE  MINV  IS  THE  STANDARD 
C  IBM  SCIENTIFIC  SUBROUTINE  MATRIX  INVERTER. 

CALL  MINV(UINV.6*D.L.M) 

C 

C  OBTAIN  THE  SOLUTION  VECTOR 
DO  69  MK»I.6 
RES(MK) »0. 

DO  69  ML-1.6 

69  RES(MK)  •A(H.)*UINVCMC»rt-3+RES<flO 


uuuu 
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EOUATION  20 

BEGIN  TO  SET  UP  THE  COEFFICIENT  MATRIX 

ua.n-i* 
uc2«  n-0. 
u(3.n-i. 

U(4>  t)*0« 

U(5.  t)*l • 

U<6,l)-0. 

ua*2)-i. 

U(I«3)*l. 

U<2.3>— 2. 

U(3.3)— 3. 

U(4#3) *2. 

U(U4)-1. 

U( 1'5) —  1 . 

U(2.5)-0. 

U(4,5>-0. 

U(S»5) “0* 

UC6.5>-0. 

UU.6>  —  1. 

U(5.6>-0. 

U  (4*  2)  ■- 1  .  '<3 --2 .  *NUH) 

U(3.2>— 3.  *0(4.2) 

UC2*2>*(3.-4.«HUM>*UC4#2> 

U(5*2)-U(3,2)>^ 

U(6#2)»U(4#2)/VF 

UC5-3)  — 3.*VF<^F 

U<6.3)-2.*vF*vF 
U(2#4)«l ,/( l.-2.*WUM> 

U(3#4)  «U(2*4) 

0(4.4)— U(2. 4) 

U(5.4)-VF«OC2.4) 

UC6.4)  — U(5*4) 

U(itlw§^4.**JF>/'<3.-2.*WF> 

U(4.6)-SFGM/'<3.-2.*HJF> 

U(3.6>— 3.*UC4.6> 

C  STORE  THE  COEFFICIENT  MATRIX  INTO  UINV. 
DO  81  1H-1-6 
DO  8t 

B1  UINV(IH-JH)-UCIH.JH) 

C  FORMULATE  THE  RIGHT  HAND  SIDE. 

A<l>*0. 

A(2)-0. 

A<3)-0. 

A(4)*0. 

A<5)-1. 

A<6)-8. 

C  INVERT  THE  MATRIX. 

CALL  MINV(UINV.6.D.L.M) 

C  OBTAIN  THE  SOLUTION  VECTOR. 

DO  89 
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RE$(rtO-0. 

DO  89  M*t.6 

RES ! MO  -A ( ft. ) *U 1 NVtrtC ,  M. )  4RES  !rtO 
GTGHC2) •!  ./!  I  ,+2.*C  l  .-NUM)/U.-2.*NUM)*VF*«£5(4)) 

C 

C  EQUATION  21 

KM GMI«<KMGMM.>*vF*2. 

GFGMl»!GFGM+l.)*VM*KMGM 

GTGMC  "GFGM*KMGM 1  *2 .  *GFGM*VmCFGM  1 

GTGMC3)  •GTGMCxCKnGm +2 .  •GFGn^VmGFGfU  ) 

C 

C  EQUATION  22 

GTGM!4)  •  l .  +VF/U . /(GFGM- 1 . )  +iKrGt»2 . ) *Vtt/<2 . *OCMGMM . ) )  ) 

C 

C  EQUATION  23 

FRAC*(  (KFGF+2. )  *VF/2 .  /!KFGM*CFGM)  *M./(1.  -GFGM) ) 

100  GTGMC5)  -GFGM+VM/FRAC 

URITECIOlRT,  425) 

425  FORMAT!//.'  TRANSVERSE  SHEAR  MODULUS'/) 

C  OUTPUT  DATA  ON  DATA  FILE 
DO  683  >1.5 

003  GTGM(J)-GTGM<J)*GM 

DO  33  J-t.5 
JJ-J+ie 

WRITE! I0URT.24)  JJ,GTGM!J> 

33  CONTINUE 

C 

C  CALCULATE  THE  MINOR  POISSON  RATIO  USING  THE  RULE  OF  MIXTURES 
C  TRANSVERSE  YOUNG'S  MODULUS  AND  EACH  OF  THE  PREVIOUSLY  CALCULATED 
C  TRANSVERSE  SHEAR  MODULI  I. 

C 

E22-ETEMU) 

URITE( 10URT.79) 

DO  75  JM.5 
JJ-J+I0 
C 

C  EQUATION  t 

P0!23-CE22/!2.*GTGM(J)))-t. 

IFCPOI23.LT.0.0.OR.POI23.GT.0.5)  GO  TO  76 
WRITE ! IOURT. 77)  JJ.P0I23 
GO  TO  75 

76  URITE<l0uRT.78>  JJ 

75  CONTINUE 

77  FORMAT!'  EOU.  I  USING  EOU.  '  .  I2.4X,  IPEU  .4) 

70  FORMAT!'  EOU.  1  USING  EQU.  ',12.'  NOT  APPLICABLE') 

79  FORMAT!//.'  MINOR  POISSONS  RATIO'/) 

C 

C  CALCULATE  THE  LONGITUDINAL  SHEAR  MODULUS 
C 

C  EQUATION  24 

GLGM!  t )  • (GFGM*!  I  ,+VF)  +VM)  /(GFGM*VMH .  +VF ) 

C 

C  EQUATION  25 

GLGM!  2 )  •  1 . /( VM+GMCFoVF) 

C 

C  EQUATION  26 

ETA ■ VF* ( GFGM- 1 . ) / ( GFGM+ 1 . ) 

GLGMC3) •! 1 ,+ETA)/! I .-ETA) 

C 

C  EQUATION  27 


CL  W>*2 .  -CFCfV  (GFGrt- 

CLGrt<4)  -GLGrt>/<2.MGFGn~ 1 .  >*VM> 


E0UAT10N  26 
FRACM./ 
CLGTK5)- 


FGrt-I.) 

i*VFy<FRftC^vn^2.) 


C  EQUATION  29 

FRAC-l.'0.-GFGm*VE/<2.*CFGM) 

GLGM(6)»GFGh+VfVFRAC 

c 

c  EQUATION  30 

ALPHA-GFGtt»(4.*VF*P I ) *P 1-4 . *VF 
ALPHA*ALPHA/<GFGrt*<PI-4.*VF>4PI*4.*^) 

GLGNC7>-.5*C(4.-PI^I*ALPHA)/4.*4.*ALPHA/'(ALPHA*(4.-P1>4PI)> 

OMTEUOURT,  325> 

325  FORMATS.'  LONGITUDINAL  SHEAR  fCDULUS'') 

C  OUTPUT  DATA  ON  DATA  FILE 
GrtG*  AflAX !  i  GF .  GM) 

GUi*AMlM|lGF,GM) 

DO  602  J-l*? 

002  GLGM(J)  •GLGM(  JMCM 

DO  34  JM#7 

IF(GLGn(J).GT.GHG.OR.GLGf1(J).LT.GLU)  GO  TO  85 
UR1TE(  IOlRT.24)  JJ.GLGM(J) 

GO  TO  34 

05  kft!TE<10t*T,63>  JJ 

34  CONTINUE 

C  MORE  DATA 

tRlTECIOCRT,  990) 

990  FORMATOV*  fDPE  DATA?#x 

1  *  INPUT  EF,EM.NUF*NUfUVF  F0RMATt5F)V 

2  •  NEGATIVE  EF  STOPS  THE  PROGRAM' ) 

READt IORED. !0)EF.EM.NUF*NUM.VF 

IF (EF .GT.0. )  GO  TO  969 

STOP 

END 
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FIGURE  CAPTIONS 


Fig.  1.  Direction  of  loading  for  the  determination  of  the  longitudinal 
Young’s  modulus  and  the  major  Poisson’s  ratio. 

Fig.  2.  vs  diagram  for  Kryo-210  superconductor* 

Fig.  3.  Normalized  plot  of  the  longitudinal  Young’s  Modulus  versus 
volume  fraction  of  the  fiber  (NbTi):  In  this  graph,  as 
well  as  the  other  normalized  comparison  graphs,  the  legend 
refers  to  the  theoretical  equations  presented  in  Appendix  I* 

Fig.  4.  e2  vs  diagram  for  Kryo-210  superconductor. 

Fig.  5.  Normalized  plot  of  the  major  Poisson’s  ratio  vs  volume 
fraction  of  the  fiber  (NbTi). 

Fig.  6.  Direction  of  loading  for  the  determination  of  the  transverse 
Young's  Modulus  and  the  minor  Poisson's  ratio. 

Fig.  7.  a2  vs  e2  diagram  for  Kryo-210  superconductor. 

Fig.  8.  Normalized  plot  of  the  transverse  Young's  Modulus  versus 
volume  fraction  of  the  fiber  (NbTi). 

Fig.  9.  vs  e2  diagram  for  Kryo-210  superconductor. 

Fig.  10.  Normalized  plot  of  vs  the  volume  fraction  of  the 
fiber  (NbTi). 

Fig.  11.  Coordinate  system  for  the  determination  of  the  longitudinal 
shear  modulus. 

Fig.  12.  Normalized  plot  of  the  longitudinal  shear  modulus  vs  volume 
fraction  of  the  fiber  (NbTi). 
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SUPERCONDUCTOR  (NbTi)  FILIMENT 


Figure  1 


IU 


Ui/UJ 


G23/G  Co 

cp.60  0.70  0.00  0.90  1.00 


08*0  0 Z.*(P 


